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Abstract 

We present several examples of quasi-exactly solvable iV-body problems in 
one, two and higher dimensions. We study various aspects of these problems 
in some detail. In particular, we show that in some of these examples the cor- 
responding polynomials form an orthogonal set and many of their properties 
are similar to those of the Bender-Dunne polynomials. We also discuss QES 
problems where the polynomials do not form an orthogonal set. 
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I. INTRODUCTION 



In last few years, the quasi-exactly solvable (QES) problems have attracted a lot of 
attention M . In particular, by now a detailed study has been made of several QES problems 
in (one particle) non-relativistic quantum mechanics in one dimension and several interesting 
features have been uncovered. However, very few QES problems in multi-dimensions or N- 
body problems in one dimension |^||] have been discussed so far. Further, to the best of 
our knowledge, no QES iV-body problem in two and higher dimensions has been discussed 
so far. The purpose of this note is to initiate a systematic study of the various aspects 
of the iV-body QES problems in multi-dimensions. In particular, one would first like to 
discover several iV-body QES problems in two and higher dimensions. Further, one would 
like to know if the corresponding polynomials form an orthogonal set or not. Besides, one 
would like to know if there is an underlying hidden algebra. In this context, it is worth 
recalling that in the case of the one body problem in one dimension, it is known that in 
most QES cases (though not all H), there is an underlying hidden algebra S7(2). Further, 
the polynomials form an orthogonal set provided the Hamiltonian can be written in terms 
of the quadratic generators of this S7(2) algebra || and that in this case one has one || or 
two sets of Bender-Dunne polynomials. 

In this paper we discuss in some detail three different iV-body QES problems namely, 
Calogero-Marchioro type iV-body problem in D-dimensions |||| in Sec. II (D > 2), iV-body 



system exhibiting novel correlations in two dimensions |10| in Sec. IIIA, and Calogero- 
Sutherland type iV-body problem in one dimension [11] in Sec. IIIB. In each case we discuss 
QES problem in case the potential is either of the sextic-type or is a sum of oscillator plus 
Coulomb type-potential. We show that whereas in the former case, the polynomials form 
an orthogonal set, in the later case they do not. We discuss the various properties of the 
orthogonal polynomials in these problems and show that they are similar to those of the 
Bender-Dunne polynomials in the one-body case. As a special case, we also discuss a self- 
dual system |T^] in which case we can obtain several more QES levels analytically. Further, 
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in this case we obtain a novel relation between the weight functions of the system. Finally, 
in Sec. IV we summarize our conclusions and point out some of the open problems. 



II. CALOGERO-MARCHIORO TYPE PROBLEM IN D-DIMENSIONS 

The N-body Hamiltonian corresponding to the Calogero- Marchioro Model in D- 
dimension for arbitrary potential V can be written as [|9] 

i=l *<i W i<j, i,j^k 1 ki 1 kj \i<j ) 

where = (r^ — r 3 -) and is the D-dimensional position vector for the i-th particle. 
On substituting the ansatz 

^=(n4J m ( 2 - 2 ) 

in the Schrodinger equation (Hip = eip) corresponding to the above system ( with h = m = 
1), it is easily shown that the equation satisfied by 0(p) is || 

9F -I- 1 

f(p) + -^-^<j>'{p) - 2V{p)4>{p) + E<j>{p) = . (2.3) 



Here E = 2e while 



1 N 

2 1 \ ~> 2 

p =^E^- 

i<3 

T D = ~ [D(JV - 1) - 2 + A D iV(iV - 1)] 



Ad = VG = \ 



\D-2f + Ag-(D-2) 
We now consider two different forms of V(p) and obtain QES solutions in each case. 



(2.4) 



A. Sextic Potential and QES problem 



Let us first consider V(p) as given by 



V(P) = - 



Bp 2 + Cp 4 + Hp 6 + — 
p 2 



(2.5) 
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In the special case when C = = H, an infinite class of exact solutions including the bosonic 
ground state have already been obtained P,|T3|. 
On substituting 



(j>{p) = Ax P (-«p 2 - Pp a )v(p) 



(2.6) 



in Eq. fl2.3|) we obtain 

v"(p) + 



2a + 2T D + 1 
P 



[Aap + 8(3p 3 ] )v'(p) + 



(E - Aa(a + T D + 1) + p 2 [4a 2 - B - 8/3(a + T D + 2)]) ??(p) = 



(2.7) 



where we have chosen, a 



c 



(3 = and F = a 2 + 2aYn . Thus a is zero or nonzero 



depending on if F is zero or nonzero. Finally, on substituting 

P n (E)p 2n 
V{P> V4"n!(n + a + T D )! 

in Eq. fl2.7|) we obtain the recursion relation satisfied by P n {E) as 



(2.8) 



P n (E) + [E - Aa (2n - 1 + a + T D )] P n -i(E) + 4(n - l)(n - 1 + a + T D ) 

4a 2 - 5 - 8/3 (2n + a + T D - 2)1 P n - 2 (E) = (2.9) 



with initial conditions P_i = and Po = 1. Using the well known theorem [14, la]], " the 
necessary and sufficient condition for a family of polynomials {P n } (with degree P n = n) to 
form an orthogonal polynomial system is that {P n } satisfy a three-term recursion relation 
of the form 



P n (E) = (A n E + B n ) P n ^(E) + C n P n . 2 {E) , n > 1 



(2.10) 



where the coefficients A n , B n and C n are independent of P, A n ^ 0, Ci = 0, C n 7^ forn > 1" 
, it then follows that {P n (E)} for this problem forms an orthogonal set of polynomials with 
respect to some weight function, uj(E). 

Let us write B in the form B = 4a 2 — 8/5(2 J + a + Td) where J is any arbitrary number. 
We shall see that when J is a positive integer, then this represents a QES system. In terms 
of J, Eq. (|2.9|) then can be rewritten as 



P n (E) + [E-Aa(2n-l + a + T D )] P n _i(£) 

-64/3(n - l)(n - 1 + a + T D ) (n - J - 1) P n _ 2 (£) = (2.11) 

from where it is clear that so long as J is positive integer, this recursion relation will reduce 
to a two term recursion relation. Thus, when J is a positive integer, we have a QES system. 
These recursion relations generate a set of orthogonal polynomials of which the first few are 

P 1 = -E + Aa(a + T D + 1) 

P 2 = E 2 -E [Sa(a + T D + 2)] + 16a 2 (a + T D + l)(a + T D + 3) 

-64/3(a + r D + l)(J-l) 

(2.12) 

It is easily seen that when J is a positive integer, exact energy eigenvalues for the first 
J levels are known. Further, when J is a positive integer, these polynomials exhibit the 
factorization property as given by 

P n+J (E) = Pj(E)Q n (E) (2.13) 

where the polynomial set Q n {E) correspond to the non-exact spectrum for this problem 
with Qo(E) = 1. For example, for J = 1, P n+ i will be factorized into P 1 and Q n and 
the corresponding QES energy level (which in this case is the ground state) is obtained by 
putting Pi = i. e. 

E l = Aa(a + Y D + 1) (2.14) 

Similarly for J = 2, P n +2 will be factorized into P 2 and Q n and the corresponding energy 
levels are, 

E l = 4a(a + T D + 2) - Ay/ a 2 + 4/3(a + T D + 1) 

E 2 = Aa(a + T D + 2) + Ay/ a 2 + 4/3 (a + T D + 1) (2.15) 

The quotient polynomials Q n {E) also form an orthogonal set as they satisfy the recursion 
relation 
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Q n (E) + [E-4a {2{n + J) - 1 + a + T D }} Q n -i(E) 

-Q4(3(n + J - l)(n + J - 1 + a + T D ) (n - 1) Q„_ 2 (£) = (2.16) 

with Qo(-E) = 1 and Q-i(E) = 0. 

The square norm of both P n and polynomials can be calculated from the recursion 
relations in Eqs ( 2.11Q and Q2.16 ) respectively. We obtain 



7^ = 64/? J] k (k + a + T D ) (J - n) 

k=l 

n 

^ = 6APl[(k + J)(k + J + a + T D ) (2.17) 

fc=i 

Note that a la the Bender-Dunne case |J the norm for the P n polynomials vanish for 
n > J while it is positive definite for n < J. This is an alternative characterization of the 
QES system. On the other hand, the norm for the Q n polynomials is always non-vanishing 
and positive definite. 

One can also calculate the weight functions uik for the set P n {E) by using the relations 
j 

Pn(E k )u k = 6 n0 , n = 0, 1, J - 1 . (2.18) 

k=l 

For example, for J =2 we obtain 

1 a 
= - + 



2 2Ja 2 + 4/3(a + T D + 1) 



uj 2 = - a (2.19) 

2 2 J a 2 + 4/3(a + T D + 1) 



Notice that to the leading order, the weight functions are inversely proportional to the 
particle number iV for a fixed D, while they are inversely proportional to the square root 
of the space dimension D, when the particle number N remains unchanged. It may also 
be noted that both the weight functions u)\ and u 2 are positive definite. Actually one can 
prove on very general grounds that for any J, the weight functions in this case will always 



be positive. As has been shown by Finkel et al. ||15|| , if the three-term recursion relation is 
of the form 
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P k+1 = (E- b k )P k - a fc P fe _! , k > (2.20) 

with a = and a n+ i =0, then (i) the weight functions are all positive (ii) the norm of the 
polynomials (7^) is positive for n < J, provided if b k is real for < k < n and a k > for 
1 < k < n. Using Eq. (|2.11 ) it is then easy to see that hence for any J, the weight functions 



and 7^ (for n < J) will be positive. 

One can also calculate the moments /i„ = / dEE n uo{E) of the weight functions. It is 
easily shown that both the odd and the even moments are non-zero. Further, to the leading 
order 

fi n = [4a(a + T D + l)] n + ■■■ . (2.21) 

Thus to the leading order, the growth rate is proportional to N 2n for fixed D and proportional 
to D" for fixed N. 

Finally, let us discuss the consequences of the anti-isospectral transformation (also termed 
as duality transformation) |T2[ in the context of our QES problem. It is easily seen that 
because of the duality transformation (x — > ix), the QES levels of the potential in Eq. 
Q2.5D are related to that of a similar potential (say V) where C (and hence a) has been 
replaced by — C. In particular, if J levels of the potential in Eq. (|2.5|) are QES levels 
with energy eigenvalues and eigenfunctions E k and ipk respectively {k = 1, 2, J), then the 
corresponding QES eigenvalues and eigenfunctions of V are given by 

E k = -E J+1 _ k , ip k {x) = ip J+1 - k (ix) (2.22) 

For example, for J =2, the QES energy eigenvalues of V are same as those given by Eq. 
( |2. 151 ) with the precise relationship being E\^% = — -E^i- The other properties of the dual 
potential are also similarly related and can be easily worked out. 



B. Self-dual Potential 



Let us now discuss the special case of C = in Eq. ( |2.5|) in which case we have a self-dual 
QES potential which has remarkably simple properties. For example, the recursion relation 
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( |2-11| ) now takes the simpler form (note that now a = C/AyH = 0) 



P n {E) + EP n ^(E) - 64/3(n - l)(n - 1 + a + T D ) (n - J - 1) P n _ 2 (S) = (2.23) 
and similarly the set Q n satisfy Eq. ( J2.16| ) with a = 0. As a result, both P n (E) and Q n (E) 



are eigenf unctions of parity. First few polynomials generated by this recursion relation are 



given by (s = 1 + a + To) 
P l = -E 

P 2 = E 2 - QA(3s(J - 1) 

p 3 = —E 3 + E [64(3s(J - 1) + 128/3(9 + 1)(J - 2)] 

P 4 = - £ 2 [64/3s(J - 1) + 128/3(s + 1)( J - 2) + 192/3(s + 2)(J - 3)] 



From the duality relation (|2.22| ) it follows that for this self-dual system, all the energy 
eigenvalues must be symmetrically distributed around E = 0. Further, in this case QES 
eigen values can be analytically calculated very simply for J = 1, 2, 5. For example, the 
eigen values for J = 3 and J = 4 are given as 



+ (64 / 3) 2 3s(s + 2)( J - 1)( J - 3) 



(2.24) 



J 



3 



8^2/3(25 + 1) 








(2.25) 



and 



J = 4 




(2.26) 
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The square norms of the polynomial sets P n and Q n s are again given by Eq. fl2.17|) . 



Since the polynomials have definite parity the odd moments vanish. The even moments of 
the weight functions in the leading order are given by 

/i 2 n~ [(l + a + T D )] n + ••■ (2.27) 

This means [i2 n is proportional to N 2n for fixed D and proportional to D n for fixed N. It is 
worth noting that in the limit (7 = (and hence a = 0), the leading term in the expansion 
of \x n in Eq. (|2.21|) vanishes, and that is why the leading behavior of /i n as function of 



(1 + a + r) when C ( and hence a) is non-zero is same as that of fj, 2n when C ( and hence 
a) = 0. 

The weight functions in this self dual QES model are related through a nice relation 

= uj+i-k (2.28) 

This can be proved explicitly using the relation in Eq. ( |2.18j ) as follows 
We consider n to be odd in Eq. (|2.18|) , i.e. 



^P 2m+1 (£ fc V fc = for k — 1, 2..., J (2.29) 

k 

Now since the P n (E) are of definite parity of E, we can write 

m 

P 2m+1 (E k ) = Y j anEl n+l , m = 0,1,2.-. (2.30) 

n=0 

when (E- independent) constant and a m — — 1. 

Putting in ( ^29|) we obtain 

J m 

J2 E a n El n+1 uj k = (2.31) 

k=l n=0 

This implies that 

E^ + W = (2.32) 
fe=i 

Now in Eq. ( |2.32| ) we use the duality relation between the eigen values of the self-dual 



system, i.e. E k = —Ej + i_ k to obtain 



E J+l-k u k = o 



(2.33) 



k=l 



Next we make a change of variable k = J + 1 — j in Eq. (|2.33|) to obtain 



(2.34) 



Subtracting Eq. ( gjp from Eq. (^32|) 



^E fc 2m+1 [^-^ J+1 _ fc ] = 

k=l 

This is true for arbitrary and hence the result, uj^ = wj+i-t- 

We now explicitly calculate the weight functions for the case of J = 3 and J 
find, 

J = 3 : 



(2.35) 



4. We 



^2 



^3 
S + l 



2(2s + i; 



2s + 2 



and for 



(2.36) 



J = 4: 



LUi = Cc>4 



1 

4 



2s + 5 



U 2 =U} 3 = - 



16s(s + 2) + 25 
2s + 5 



16s(s + 2) + 25_ 

Not surprisingly, the weight function satisfy relation ( |2.28|) and X^^i^fc — 1- 



^2.37) 



C. An unusual QES Problem 



Let us now discuss an unusual QES problem. To that end consider the following potential 

V(p) = £[SY-£ + $] (2-38) 
z p p 
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It is worth noting that for C = as well as for B = 0, a class of energy levels can be 
analytically obtained [§,0- For the special case of N = 2 i.e. 2-anyons in a magnetic field 
and experiencing Coulomb potential . This problem has been discussed by Myrheim et al. 
fTTfl . We want to show that when both of them are non-zero, it is an example of an unusual 
QES system in the sense that for one parameter family of potentials (when there is a specific 
relation between B, C and F), one always obtains one QES eigenvalue. 
On substituting 



<f>(p) = p a ex V (-Bp 2 /2) V {p) 



(2.39) 



in Eq. ( |2.3|) we obtain 

"2a + 2T D + 1 



V"(P) + 



P 



(2Bp) 



v'(p) + 



E-2(a + T D + 1)B 



C 



■nip) = o 



(2.40) 



where we have chosen, F = a 2 + 2aTo . Next we substitute 



V(p)=E P n(E) P n 

n 

111 Eq. ( |2~40D and obt ain the recursion relation satisfied by P n {E) as 



(2.41) 



n{2a + 2T D + n)P n {E) + CP n ^(E) + [E - 2(n - 1 + a + T D )]P n ^ 2 {E) = (2.42) 

with initial conditions P_i = and Pq = 1. 

From the recursion relation it follows that the QES levels are obtained in case the coeffi- 
cient of P n -2(E) vanishes, i.e. the energy of the QES levels is given by E = 2{a + Yn +n+ 1) 
where n = 1,2,.... Further, in that case B, C and F are not independent. For exam- 
ple, for n = 1, the relation is C 2 = 2(2a + 2r^ + 1)B while for n = 2 the relation is 
C 2 = 4(4a + 4Td + 3)B. The QES level corresponds to an excited (ground) state depending 
on if C > (<)0. On comparing Eq. ( |2.42j ) with Eq. ( [2.10 ), we conclude that the polynomials 
in this case do not correspond to orthogonal polynomials. 
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III. OTHER N-BODY PROBLEMS 



In this section we briefly discuss two other many body problems and we will see that in 
both of these cases we have two QES problems which are very similar to those discussed in 
the last section. 



A. Novel Correlations for N- particles 

Recently, a class of exact solutions including the bosonic ground state has been obtained 
for a many body Hamiltonian in two dimensions such that the wave-functions have a novel 
correlation (of the form Xy = XiUj — Xjt/i) built into them [Hj. Here we would like to show 
that QES states with novel correlations can also be obtained in this case. To that purpose 
we start from the Hamiltonian 



H 



fc2 i-2 r 2 i-2 

— — W 2 + — — n V + — — n V j ' k 

2m ^ * 2m yi 4? X% 2rn ^ X u ■ X ik 



L K i i i 



Yr 2 



On substituting the ansatz 



A' 



ij = Y[\x ij \^( P ) 



P 



E-Bp 2 - Cp A - Hp 6 - 







Where 



(3.1) 



(3.2) 



in the Schrodinger equation Hip = eip, (m = H = 1) one can show that satisfies the 
equation 



(3.3) 



2A + 1 = 2N - 1 + 2gN(N - 1) 



E = 2e; g 1 = g{g - 1) , g 2 = g 2 



(3.4) 



On substituting 
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<j)(p) = p a exp {-ap 2 - pp 4 )v(p) 



(3.5) 



in Eq. ( |3.3|) we obtain 



v"(p) + 



2a + 2 A + 1 
P 



- (Aap + 8(3p 3 



v'(p) + 



E - Aa{a + A + 1) + p 2 (4a 2 -B- 8(3{a + A + 2))] r]{p) = 



(3.6) 



where we have chosen, a 



c. p _ Vh and F _ fl 2 + 2a A . Now notice that this 
equation is similar to the corresponding equation in the last section and hence the rest of 
the analysis of that section will go through step by step in this case. Thus this is another 
example of QES iV-body problem in two dimensions where the polynomials satisfy most 
of the properties of the Bender-Dunne polynomials. Similarly, one can also consider the 
admixture of the oscillator and Coulomb type potential and as in the last section, again 
obtain similar QES solutions. In particular, in this case the polynomials do not form an 
orthogonal set. 



B. N-body QES problem in 1-dimension 

Finally, we would like to briefly discuss the Calogero-Sutherland type iV-body QES 
problems in one dimension ||,|3|]. 

The Hamiltonian corresponding to the N-body QES problem in 1-dimension is given by 

tfE*. 2 l 3 + ^? (3.7) 
We substitute the standard ansatz for the wave function ip, 

iP = Z x+ ^(j)(p) (3.8) 

where Z = rii<j( s i — x j), p 2 — J2%h ^ = |\A + ^9 i n = e ^P = m = 1) an d obtain 
the ^-equation 
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2A + 1 
P 



where E = 2e and 



+ 



E-Bp 2 - Cp 4 - Hp 6 



F 



(3.9) 



2A + 1 = A-l + A(iV-l)A 



(3.10) 



This equation is very similar to the Eq. (|3.3|) and hence the rest of the analysis goes through 
as in the last section, i.e. the polynomials of the QES problem indeed form an orthogonal 
set and share most of the properties of the Bender-Dunne polynomials. Similarly, we can 
show that even in this case the other QES solution corresponding to the admixture of the 
oscillator and Coulomb potential also exists, and in this case, the polynomials do not form 
an orthogonal set. 



IV. SUMMARY AND OPEN PROBLEMS 

In this paper, we have obtained QES solutions of a number of TV-body problems in 
one,two and higher dimensions. We have shown that with sextic type potential, in all these 
cases one has a set of orthogonal polynomials which satisfy almost all the properties as 
satisfied by the Bender-Dunne polynomials. However, the hidden algebra is not clear in these 
cases. In this context, recall that for the case of one particle in one dimension experiencing 
the sextic potential, the underlying symmetry algebra is Sl(2). In particular, in that case 
one can express the Hamiltonian in terms of the quadratic (and linear) generators of Sl(2). 
Further, it will be nice to discover some other QES many-body problems. Finally, it would 
be most interesting if one can find some application of these QES A-body problems. In this 
context it is worth noting that one application has already been found by one of the present 



author (AK) and Jatkar [0. In particular, they have considered the model discussed in 
Sec. IIIB (with F =0) and have shown that the square of the ground state wave function of 
this model is related to the matrix model corresponding to branched polymers. It will really 
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be interesting, if one can find some other application, of the N-body QES problems in two 
and higher dimensions. 
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